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Preface 


My  Initial  motivation  for  this  thesis  was  to  develop  and  evaluate 
a  new  design  procedure.  I  would  be  less  than  honest  If  I  failed  to  also 
admit  that  I  was  further  motivated  to  complete  a  thesis  which  would  meet 
the 'requirements  for  graduation.  Somewhere  along  the  way  motivation 
shifted.  I  hope  the  many  long  hours  of  work  sumnarlzed  In  this  thesis 
might  be  of  use  to  others  working  on  the  landing  problem.  It  Is  then, 
to  this  end,  that  I  have  written  this  thesis.  I  have  Included  the 
equations  and  the  numerical  values  which  I  used  so  others  will  be  able 
to  reproduce  my  results.  I  have  also  tried  to  write  this  thesis  In  a 
clear,  understandable  manner  for  the  benefit  of  the  prospective  reader. 

Those  who  have  written  a  thesis  can  appreciate  the  many  long  hours 
of  work  required  and  the  fact  that  a  thesis  Is  not  just  an  Individual 
effort.  I  would  like  to  express  my  sincere  appreciation  to  those  who 
helped  make  this  thesis  possible.  I  especially  want  to  thank  my  advisor. 
Major  James  0.  011  low,  for  his  many  hours  of  help  and  encouragement.  His 
cheerful  Npress-onM  attitude  turned  a  dreaded  task  Into  a  learning  exper> 
lence.  Major  Dlllow  deserves  a  majority  of  the  credit  for  this  study. 

I  would  also  like  to  thank  my  sponsor,  Robert  Huber,  and  Ronald  Anderson, 
both  of  the  Air  Force  Flight  Dynamics  Laboratory,  for  their  help. 

Finally,  I  wish  to  express  my  deep  appreciation  to  my  wife,  Joyce, 
for  her  patience,  understanding,  and  encouragement,  and  for  typing  the 
rough  draft  of  this  thesis.  Her  sacrifices  during  not  only  this  thesis 
but  the  last  four  years  have  helped  Immeasurably  In  qy  education.  With* 
out  her  encouragement,  I  would  not  have  come  this  far. 

Jerry  0.  Pflecger 
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LANDING  APPROACH  AUTOMATIC  FLIGHT  CONTROL  SYSTEM  DESIGN 
VIA  REDUCED  ORDER  OPTIMAL  CONTROL  LAW 

I.  Introduction 


Background 

A  new  National  Landing  System  is  being  developed  to  provide 
improved  landing  guidance  information  to  aircraft.  One  approach  which 
has  shown  promise  involves  the  use  of  a  microwave  scanning  beam  to 
provide  landing  aircraft  with  glide  path  deviation  and  azimuth  infor¬ 
mation.  James  D.  Dillow  conducted  a  study  (Ref  1)  to  analyze  the  pitch 
plane  data  rate  requirements  for  the  proposed  microwave  scanning  beam 
system.  The  data  rate  study  used  an  optimal  model  for  the  aircraft 
flight  control  system. 

As  a  part  of  the  data  rate  analysis,  a  digital  computer  program 
was  developed  to  implement  the  landing  approach  model  and  to  automatically 
compute  the  probability  of  a  missed  approach  as  a  function  of  the  data 
rate  of  the  glide  slope  deviation  information.  The  computer  program 
takes  as  inputs  the  data  rate;  the  aircraft  stability  derivatives; 
nominal  longitudinal  airspeed;  glide  path  angle;  atmospheric  disturbances 
(gust,  headwind,  windshear);  guidance  noise  parameters;  constraints  on 
the  control  activity;  and  the  tolerances  on  the  aircraft  variables  that 
define  a  missed  approach  (i.e.,  the  "window  dimensions").  The  program  is 
structured  o  that  changes  In  the  system  equations  can  be  easily  Implemented 
in  the  program.  Changes  In  the  system  equations  result  from 

1.  Changes  in  the  aircraft  equations  of  motion  (for  example, 
Inclusion  of  flexure  modes  and  effect  of  sensor  location). 
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2.  Addition  or  deletion  of  control  points  (i.e.,  autothrottle, 
direct  lift  control,  etc.). 

3.  Changes  in  the  measurement  model  or  sensor  complement  (for 
example,  consideration  of  sensed  normal  acceleration,  sensed 
glideslope  deviation  rate,  sensed  longitudinal  airspeed,  and 
noise  in  the  continuous  measurements). 

The  results  obtained  by  Dill ow  (Ref  1)  using  the  optimal  model  were 
compared  with  results  obtained  using  conventional  analysis  techniques 
(Ref  2).  This  comparison  validated  the  optimal  model. 

One  of  the  "by-products"  of  the  digital  computer  program  used  in 
the  data  rate  analysis  is  a  "full-blown"  optimal  control  law  with  feedback 
gains  and  Kalman  Filter  gains.  This  control  law  Is  optimal  In  the  sense 
of  approach  performance  and  accounts  for  control  authorities  that  would 
be  Imposed  on  the  automatic  flight  control  system.  The  purpose  or  object 
of  this  study  was  to  determine  If  the  optimal  model  and  the  optimal  control 
law  from  the  data  rate  analysis  could  be  used  as  a  design  tool.  The  study 
was  prompted  by  the  success  of  the  data  rate  analysis  and  the  fact  that 
the  "optimal"  was  optimal  with  respect  to  a  meaningful  measure  of  per¬ 
formance  while  accounting  for  limitations  in  control  activity.  The 
digital  computer  program  used  in  the  data  rate  analysis  offers  several 
advantages  over  classical  techniques  in  control  design.  The  program 
structure  makes  changes  in  the  system  equations  a  matter  of  changing  data 
cards  rather  than  redrawing  root  locus  or  Bode  plots  as  is  done  with 
classical  techniques.  There  is  no  need  to  go  through  tedious  loop  closure 
procedures.  The  time  required  to  evaluate  the  effect  of  one  or  more 
changes  is  much  less  than  that  needed  when  using  classical  techniques. 
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ftilti-input,  mjlti -output  control  systems  are  also  handled  Much  easier 
using  the  data  rate  analysis  computer  program. 

Hypothesis 

The  optimal  feedback  matrix  F*,  from  the  data  rate  analysis, 
can  be  used  in  a  systematic  procedure  to  develop  a  practical  control 
law.  The  systematic  procedure  is  subject  to  one  constraint.  When  the 
data  rate  analysis  program  develops  the  optimal  control  law,  certain 
variables  such  as  pitch,  pitch  rate,  normal  acceleration,  and  nominal 
longitudinal  airspeed  are  considered  to  be  measured  "on  board"  the 
aircraft.  The  data  rate  analysis  program  further  assumes  that  glide 
path  deviation  is  measured  on  the  ground  and  transmitted  to  the  aircraft. 
Therefore,  flight  path  deviation  information  is  in  sampled  data  form  for 
the  data  rate  analysis.  The  constraint  is  that  the  continuously  measured 
variables  and  sampled  data  measurements  used  in  the  data  rate  analysis 
mst  also  be  used  as  feedbacks  in  the  practical  control  law.  Any  control 
law  not  employing  all  these  feedback  gains  will  be  unstable  or  display  poor 
flying  c r  poor  ride  qualities.  These  poor  qualities  show  up  as  large 
pitch  attitude  deviations,  high  control  actuator  rates,  high  probability 
of  missed  approach  or  slow  settling  time  when  the  aircraft  Is  perturbed. 

Investigation 

The  object  of  this  Investigation  is  to  determine  a  design  procedure 
starting  with  the  "optimal*'  feedback  gains  of  the  data  rate  analysis  and 
ending  with  a  practical  (and  realizable)  automatic  flight  control  system 
for  an  aircraft  performing  the  landing  approach  task  and  to  compare  the 
resulting  suboptimal  automatic  flight  control  system  against  some  "standard." 
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The  measures  of  “goodness"  for  the  evaluation  are  the  probab111t>  of 
missed  approach,  PMA  (minimized);  the  rms  control  activity  (within  bounds); 
and  rms  pitch  attitude  o0  (ride  quality).  System  response  to  a  pitch 
rate  Impulse  will  also  be  checked  for  low  damping.  If  the  damping  is  too 
low,  the  system  will  be  slow  In  damping  out  oscillations.  These  oscilla¬ 
tions  are  not  evident  in  the  rms  performance  measures. 

The  "standard'  used  as  a  basis  for  comparison  was  designed  by 
Systems  Technology,  Inc.  (Ref  2).  It  represents  an  advanced  automatic 
flight  control  system  designed  with  classical  control  techniques.  The 
"standard”  is  not  currently  in  use,  but  rather  represents  a  system  of 
the  quality  deemed  necessary  to  meet  the  Category  II  landing  requirements. 

In  order  to  compare  results  with  the  "standard"  a  common  basis 
must  be  established.  The  following  ground  rules  have  been  established 
to  arrive  at  that  basis: 

1.  The  glldeslope  deviation  d  Is  assumed  to  be  measured  contin¬ 
uously  in  the  "standard"  and  sampled  at  a  rate  of  6  samples/sec 
In  the  suboptlmal  control  system. 

2.  The  sample  data  glldeslope  deviation  measurement  Is  assumed 
to  contain  noise  where  the  noise,  which  Is  superimposed  on 
the  true  glldeslope  deviation,  accounts  for  the  effects  of 
fluctuation  noise  (due  to  sampled  data  measurement)  and  white 
noise.  This  accounts  for  using  a  scanning  beam  landing  guidance 
system.  The  noisy  measurement  is  filtered  with  a  low  pass 
filter  which  Is  described  in  detail  later.  When  airspeed  Is 
used.  It  Is  assumed  that  any  noise  which  might  be  on  the 
measurement  Is  filtered  out  by  the  measurement  system.  The 
measurement  system  for  airspeed  is  approximated  by  a  first- 
order  lag. 
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3.  The  atmospheric  environment  is  one  of  severe  turbulence. 

4.  The  microwave  scanning  beam  has  low  beam  noise. 

5.  The  height  Ad  of  the  landing  "window"  is  ±12  ft. 

6.  Steady  winds  and  wind  shear  are  both  assumed  to  be  zero. 

7.  The  measured  variables  are  pitch,  pitch  rate,  normal  accelera¬ 
tion.  and,  in  special  cases,  longitudinal  airspeed. 

The  Individual  ground  rules  are  explained  in  greater  detail  as  they 
are  encountered  later  In  this  report.  The  procedure  is  explained  here 
In  general  and  In  detail  in  Chapters  V  and  VI. 

A  basic  requirement  for  the  procedure  Is  that  it  be  a  "clean"  and 
logical  process.  The  starting  point  Is  the  full  optimal  control  law  using 
the  optimal  feedback  matrix  F*  from  the  data  rate  analysis.  The  goal  Is 
a  reduced  order  control  law  employing  feedback  of  the  measured  variables 
(states)  only.  Feedbacks  on  the  measured  states  are  referred  to  In  this 
report  as  desired  feedbacks.  Since  the  optimal  feedback  F*  includes  some 
feedback  gains  for  variables  not  measured  (1 .e . .  wind  gust,  elevator 
position,  etc.),  the  procedure  must  systematically  remove  the  feedbacks 
on  the  unmeasured  states.  The  unmeasured  state  feedbacks  are  referred  to 
as  undesired  feedbacks.  Once  the  undesired  feedbacks  have  been  removed 
and  only  the  desired  feedbacks  remain,  the  logical  approach  would  be  to 
try  to  further  simplify  the  control  law  by  reducing  the  number  of  remaining 
feedbacks.  In  the  hypothesis,  it  was  stated  that  removal  of  any  desired 
feedback  would  produce  unacceptable  results.  The  next  step  in  the  pro¬ 
cedure  is  to  validate  this  constraint.  The  final  step  involves  placing 
noise  on  those  measurements  (glldeslope  deviation)  assumed  to  have  addi¬ 
tive  noise  present.  The  measurement  and  the  noise  are  then  filtered 
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through  a  low  pass  filter  and  the  filtered  measurement  used  for  the 
practical  control  law.  The  practical  control  law  is  developed  without 
"fiddling"  with  feedback  gains  or  use  of  compensators. 
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1 I .  Mathematical  Model 


In  this  section  the  mathematical  model  of  the  landing  approach 
process  is  described.  This  model  Is  used  to  develop  the  automatic  flight 
control  system.  The  assumptions  used  in  developing  the  model  are 
described  and  the  equations  used  In  this  study  are  given.  The  develop* 
ment  of  the  equations  Is  found  In  Ref  1.  The  system  of  equations  Is  a 
set  of  stochastic  differential  equations  and  accounts  for  the  aircraft 
dynamics,  the  atmospheric  environment,  the  landing  guidance  system  data 
rate  and  errors,  the  aircraft  onboard  sensors,  and  the  flight  control 
system  capabilities.  The  approach  performance  Is  determined  by  the 
probability  of  missed  approach. 

The  Aircraft  Equations  of  Motion 

The  aircraft  Is  assumed  to  be  adequately  represented  by  a  set 
of  perturbed,  linear  differential  equations  of  motion.  This  Is  because 
the  landing  approach  task  basically  requires  tracking  a  fixed  rectilinear 
path  In  space  with  disturbances  induced  by  the  atmospheric  environment 
(gusts)  and  measurement  errors  Induced  by  sensor  noise  and  noisy  or 
erroneous  guidance  measurements.  It  is  further  assumed  that  the  longi¬ 
tudinal  and  lateral  equations  of  motion  are  uncoupled.  This  assumption 
Is  justified  for  trimmed  flight  with  small  perturbations  in  Ref  3.  It 
Is  assumed  that  the  longitudinal  motion  variables  dominate  the  prob¬ 
ability  of  missed  approach  (Ref  1),  thus  only  these  equations  are 
considered.  The  variables  considered  In  the  longitudinal  equations  of 
motion  are  pitch  attitude  e,  pitch  rate  q,  longitudinal  perturbed 
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velocity  u;  normal  perturbed  velocity  w;  and  glides  lope  deviation  d. 
Additional  staves  are  Introduced  as  needed  to  account  for  lags  between 
cc—anded  control  Input  and  the  resulting  force  or  notion  generated, 
additive  finite  bandwidth  noise,  or  sensor  lags. 

The  specific  equations  presented  here  represent  the  system  nodel 
for  the  basic  DC-6  aircraft  and  are  adopted  fraa  Ref  1.  They  are* 

u  «  V  ♦  X*w  -  ge  cos  y,  -  Xyiig  -  X*wg  +  X6e«e  +  ^^th 
w  *  Z||U  ♦  Z^w  ♦  Utq  -  g0  sin  y#  ~  ZyUg  -  Z*wg  +  Z^fie  + 
q  -  Hyu  +  Myw  +  Mgq  -  HyUg  -  MyWg  +  n^6e  * 

The  elevator  and  thrust  responses  are  nodeled  by  first-order  lags. 
The  respective  differential  equations  are 

«e  *  “  6e  +  7"  **c 
■e  'e 

4«i  •  -  ^  «th +  ^  «t»t 

where  6e{.  and  6th  are  commanded  inputs,  6e  and  6^  are  the  resulting 
c 

elevator  position  and  throttle  setting,  respectively,  and  the  values  for 
Te  and  T*  are  taken  to  be  .06666  sec  and  1.0  sec,  respectively. 

Gusts 

The  gust  disturbance,  taken  from  Ref  1,  Is  assumed  to  have  two 
Independent  components.  These  c  wponents  are  the  longitudinal  gust 


*Z*  and  My  were  taken  to  be  zero  In  these  equations. 
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velocity  Ug  and  the  normal  gust  velocity  Wg.  These  are  described  by  the 
following  first-order  stochastic  differential  equations: 

“g  ■  ‘  +  5ug 

Wg  -  -  WWgWg  +  5Wg 

where  uug  and  are  the  half  power  or  break  frequencies  for  the  longi¬ 
tudinal  and  normal  gusts,  respectively;  Zu  and  are  zero  mean, 

9  9 

Gaussian  amplitude,  white  noise  processes.  The  statistics  of  £Ug  and 
£w  are  given  by 

wg 

E(Cug(t)5Ug(s)}  ■  a^gaigitt-s) 

E<{Wg(t)CWg<S)>  - 

E<eUg<t)e„g(s)j  ■  o 

where  ou  Is  the  rms  longitudinal  gust  Intensity  and  ow  Is  the  rms 
normal  gust  velocity;  6{t-s)  Is  an  Impulse  function. 

The  parameters  used  to  describe  the  gust  disturbance  are  given  In 
Table  I.  These  values  were  taken  from  Refs  1  and  2.  The  rms  gust  Inten¬ 
sities  used  here  represent  severe  turbulence.  Severe  turbulence  was 
considered  because  previous  studies  (Refs  1  and  2)  Indicated  the  landing 
approach  task  required  landing  guidance  and  flight  control  systems  to 
suppress  the  effects  of  gusts. 


9 


TABLE  I 

Atmospheric  Disturbance  Parameters— DC-8 


0U  *  10.000  ft/ sec 
9 

Mu  -  0.340  rad/sec 
9 

a w  *6.500  ft/sec 
w9 

Mw  *  3.950  rad/sec 
9 


Measurements 

The  measurement  model  was  developed  In  Ref  1  and  accounts  for 
the  two  different  types  of  measurements  associated  with  a  low  visibility 
landing  approach  with  a  scanning  beam  guidance  system.  The  first  part  of 
the  measurement  model  accounts  for  the  glldeslope  deviation  measurements 
derived  from  the  scanning  beam  system,  and  as  such,  are  considered  to  be 
sample  data  measurements.  These  measurements  are  usually  associated  with 
guidance  measurements  or  guidance  Inputs  to  the  flight  control  system. 

The  other  measurements  considered  In  the  measurement  model  are  those 
usually  associated  directly  as  a  part  of  the  flight  control  system  and 
are  made  continuously  on  board  the  aircraft.  These  measurements  may 
Include  sensed  pitch  attitude,  pitch  rate,  airspeed,  normal  acceleration,  etc. 

The  model  for  the  glldeslope  deviation  measurement  as  derived  from 
the  scanning  beam  guidance  system  Is  taken  from  Ref  1.  This  model  for 
measured  glldeslope  deviation  accounts  for  three  error  components  In  the 
difference  between  the  measured  value  of  glide  slope  and  the  true  reference 
landing  glide  slope.  These  components  include  a  fixed  bias,  fluctuation 
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noise  associated  with  the  sampled  data  feature  of  the  measurement,  and 
white  noise  representing  the  very  broad  band  thermal  noise  In  the  circuits 
and  random  electromagnetic  noise. 

The  fixed  bias  Is  taken  as  a  zero  mean,  Gaussian  random  variable. 

It  is  not  directly  Included  as  a  part  of  the  measurement  since  It  Is 
assumed  that  this  error  cannot  be  detected  without  some  other  external 
reference.  Thus  the  true  glides  lope  plus  the  fixed  bias  Is  taken  as  the 
reference  glldeslope  track  for  the  aircraft.  Its  effect  Is  accounted  for 
by  computing  a  root-sum-squared  glldeslope  deviation  using  the  fixed  bias 
and  the  rms  glldeslope  deviation  due  to  all  other  disturbances  and  measured 
errors. 

The  fluctuation  noise  y2  Is  modeled  by  a  first-order,  Gaussian, 
shaped  noise  of  the  form 

y2  -  -  ufny2  +  £fn 

where  u>fn  Is  the  half-power  frequency  of  the  fluctuation  noise  and  £fn 
Is  a  zero  mean,  Gaussian  amplitude,  white  noise  process.  Furthermore, 

EUfn(t)£fn(s)}  .  2ufnOfn5(t-s) 

where  Ofn  Is  the  rms  fluctuation  noise.  The  break  frequency  Is  taken 
to  be 

_  2.8  . 

wfn  -  —  rad 

U 

where  To  Is  the  Information  update  Interval  (or  1/To  Is  the  data  rate). 

Thermal  noise  and  random  electromagnetic  noise  are  modeled  by  zero 
mean,  Gaussian  amplitude,  white  noise  n.  It  Is  assumed  that  n  is  statis¬ 
tically  Independent  of  y2. 
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The  coDlntd  neasurtd  glldoslopo  deviation  from  the  true  reference 
glldoslopo  Is  represented  pictorial ly  In  fig.  1.  The  measured  deviation 
from  dm  reference  glldoslopo  track  y  Is  given  by 

y  ■  d  ♦  ya  ♦  n 

The  measured  glldoslopo  deviation  y  Is  sampled  at  the  Information  wpdau 
Intervals  and  y(nT«)  Is  the  naasurod  sample  data  glldoslopo  deviation  on 
the  tine  Interval  nT«  it  <  (n*1)T«,  Mhere  T«  Is  the  sampling  Interval  and 
1/Tt  Is  the  sanpled  data  rate. 


s 

s 

s 


Fig.  1.  (alldeslope  Deviation  Measurenent  Model 

In  addition  to  the  sampled  data  measurement  of  glldeslope  deviation, 
it  is  assumed  that  certain  aircraft  variables  are  continuously  measured. 

The  mathematical  model  of  Ref  1  was  developed  to  Include  consideration 
for  continuously  sensed  motion  variables;  however,  the  following 
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restriction  ms  laposed  on  the  continuous  aeasureaents  considered  In 
Ref  1: 

•  Only  those  variables  appearing  explicitly  «s  stetes  In  e 
first-order  state-variable  repretentetton  of  the  systea 
equations,  x  ■  Ax  ♦  Bu  +  C* »  cm  be  aeasured  continuously. 
Furthermore,  ttiose  aeasureaents  ere  perfect,  1  .e. ,  not 
noisy. 

The  reeson  for  this  restriction  Is  discussed  In  Ref  1.  This  restric¬ 
tion  dots  not  generally  Halt  the  type  of  continuous  measurements  th«t 
cen  be  considered  but  In  soae  esses  It  requires  "fiddling  *  round*  with 
the  stele  equetlons  to  get  then  In  e  suitable  fora  so  that  the  restric¬ 
tion  Is  satisfied;  for  example,  when  It  was  desired  to  use  airspeed 
feedback  In  the  control  law.  In  order  to  consider  longitudinal  airspeed 
(longitudinal  velocity  with  respect  to  air  aass)  as  a  sensed  variable, 
it  is  convenient  to  Introduce  the  longitudinal  airspeed  uM  as  a  state 
in  the  aircraft  equations.  Let  Ug  denote  the  Inertial  velocity  of  the 
air  aass.  Then 

was  *  u  *  ug 

and  the  tlae  derivative  of  the  longitudinal  airspeed  Is 

’  «  *  ig 

The  aircraft  equations  of  aotlon,  after  replacing  longitudinal  velocity 
with  longitudinal  airspeed,  are 

“as  "  Vas  *  V  ’  96  co*  Y#  *  ^“g  *  Vg  ♦  *«f«e  *  *«th«th 
*  "  Zuuas  *****  u#<)  ‘  9°  *1n  Y*  ’  tyg  *  *  Z«th4th 
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« •  V..  ♦  V  ♦  V  -  Vfl  *  V* 4  "«th6th 


In  «  slallar  sinner,  noisy  aeasureaents  of  aircraft  sotlon  variables 
can  bo  considered.  This  Is  dona  by  Introducing  tho  appropriate  linear 
stochastic  differential  equations  describing  the  aeasureaent  noise  and 
Introducing  a  new  state  equation  representing  the  sue  of  the  sensed  sot  Ion 
variable  and  the  aeasureaent  noise. 


mrziiiinriiinmrciir 


The  equations  of  notion  for  the  aircraft  along  with  the  control 
lags,  gust  equations,  and  neasurenent  equations  sake  up  a  set  of  dif¬ 
ferential  equations  ehlch  describe  the  systea.  The  systea  differential 
equation  Is 

x(t)  •  A*(t)  ♦  Bu(t)  ♦  C*(0  (1) 
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Tha  *  Utility  darlvatlvas  and  aircraft  paraaatars  ara  Ukan  froa  Raf  1 
and  ara  Includad  h ara  In  Tabla  II.  Tha  systaa  block  dlagraa  It  shown 
In  Fig.  2. 

Haasuras  of  ’Coodnus* 

Throa  aaasuras  of  *goodnass"  wara  usad  to  avaluata  a  posslbla 
flight  control  systaa.  Tha  first  Is  tha  probability  of  alssad  approach 
which  Is  usad  as  tha  prlairy  aaasura  of  approach  parforaanca.  Tha 
sacond  Is  tha  ras  control  ctlvlty*  and  tha  third  Is  tha  ms  pitch  angla. 

Tha  aaasura  of  approach  parforaanca  usad  In  this  study  It  basad 
on  tha  assuaptlon  that  If  cartaln  praspaclflad  aircraft  varlablas  ara 
within  a  glvan  tolaranca  at  tha  daclslon  altltuda,  than  tha  landing 
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TABLE  II 

DC-8  Stability  Derivatives  and  Aircraft  Parameters— Landing  Approach 


Xu. 

-  0.03730  1/sec 

X*- 

0.13600  1/sec 

Zu- 

-  0.28300  1/sec 

-  0.75000  1/sec 

Mu  " 

0.0  1/sec-ft 

M«- 

•  0.00461  1/sec-ft 

• 

-  0.59400  1/sec 

V 

0.0  ft/rad-sec* 

S»* 

0.10600  ft/percent  r pm -sec 

V 

-  9.25000  ft/rad-sec* 

‘‘w,* 

-  0.00097  ft/percent  rpm-sec 

V 

•  0.92300  1/rad-sec* 

%,* 

0.00007  1/percent  rpm-sec* 

u.  • 

228  ft/sec 

Y*  • 

•  3  deg 

Te* 

0.06666  sec 

1.0  sec 
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Block  Dlagraa  of  the  Landing  Approach  Model 


proceeds  to  touchdown.  Otherwise  •  missed  approach  has  occurroo  and  a 
"go  around"  Is  executed.  The  tolerances  used  to  specify  a  successful 
approach  are  called  the  window.  With  this  definition  of  •  missed  approach, 
the  probability  of  a  missed  approach  (PMA)  Is  a  quantitative  and  computable 
measure  of  landing  performanca.  The  definition,  a  complete  description, 
and  a  method  for  computing  PMA  are  found  In  Ref  4. 

The  current  FAA  Category  II  landing  accuracy  criteria  of  i72  ft 
lateral  deviation  and  ±12  ft  vertical  deviation  are  used  to  deflno  the 
window.  The  window  Is  shown  graphically  In  Fig.  3.  Since  the  results 
of  Ref  2  Indicate  the  vertical  errors  dominate  the  probability  of  missed 
approach  (with  this  window  definition),  only  the  vertical  window  dimen¬ 
sion  was  used  In  this  study  to  define  a  missed  approach. 


Fig.  3.  The  Recommended  Category  II  window 

After  minimizing  the  probability  of  missed  approach,  the  rm 

elevator  actuator  rate  end  the  ram  pitch  angle  are  checked  to  assure 
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that  both  fall  within  acceptable  Holts.  The  rms  elevator  actuator  rate 
serves  as  a  Measure  of  the  control  authority  required  by  the  automatic 
flight  control  system.  Limits  are  generally  set  on  the  control  authority 
of  an  automatic  flight  control  system  with  a  margin  of  safety  to  allow 
sufficient  manual  authority  to  override  a  hard-over  failure  of  the  auto¬ 
matic  system.  The  control  constraint  used  as  an  upper  bound  on  activity 
was  .15  rad/sec  for  the  rms  elevator  actuator  rate.  This  value  was  never 
exceeded  by  any  of  the  cases  considered.  By  relaxing  or  tightening  the 
ms  actuator  rate  constraints ,  possible  trade-offs  can  be  evaluated  In 
terms  of  the  change  in  the  probability  of  missed  approach  for  a  given 
change  In  the  rms  actuator  rate. 

The  ms  pitch  angle  was  considered  as  a  measure  of  the  "ride" 
qualities  associated  with  the  automatic  flight  control  system  and  hence 
was  used  to  evaluate  each  automatic  flight  control  system  configuration. 

The  upper  limit  was  set  at  6*  (Ref  1).  This  figure  was  primarily  used 
to  detect  an  automatic  flight  control  system  design  which  results  in  a 
high-frequency  oscillation  In  glldeslope  deviation.  There  are  cases 
where  an  automatic  flight  control  system  could  have  a  low  probability  of 
misted  approach  and  reasonable  actuator  rate,  but  the  oscillatory  tendencies 
would  show  up  by  a  marked  Increase  In  the  rms  pitch  angle.  The  probability 
of  missed  approach  and  rms  elevator  actuator  rate  appears  reasonable  in 
these  casts  because  oscillations  around  the  mean  tend  to  cancel  out  In 
the  ms  measures. 


20 


.  •  •  %»*  v  sWs-«< 


.Vi 


III.  Optimal  Feedbacks 


The  technique  used  to  develop  the  optimal  feedback  gains  Is  pre- 
sented  in  this  chapter  for  completeness  and  reader  convenience.  This 
coverage  Is  not  Intended  to  be  mathematically  rigorous.  A  more  In-depth 
treatment  Is  found  In  Refs  1  and  4. 

In  developing  the  model  for  the  automatic  flight  control  system, 
steady  winds  and  wind  shear  were  taken  to  be  zero.  The  steady  wind  is  a 
time  Invariant  wind  assumed  to  be  blowing  horizontally  with  the  ground 
and  may  be  either  a  head  wind  or  tall  wind.  With  respect  to  the  aircraft. 
It  has  a  steady  longitudinal  and  normal  component  depending  on  the  air¬ 
craft  attitude  with  respect  to  the  horizon.  The  wind  shears  account  for 
a  gradient  In  the  wind  Intensity  as  a  function  of  altitude.  The  wind 
shear  wsh  Is  represented  by  the  equation 


"sh 


,° 

^UOO-h), 


hi  200  ft 

100  ftihi200  ft 


where  sh  denotes  the  linear  rate  at  which  the  wind  velocity  changes  as 
a  function  of  altitude;  sn  may  be  positive  or  negative.  With  respect  to 
the  aircraft,  the  wind  shear  has  a  longitudinal  and  normal  component 
depending  upon  the  aircraft  attitude  with  respect  to  the  horizon. 

The  control  law  representing  the  automatic  flight  control  system 
Is  developed  for  the  linear  Gaussian  system  of  equations. 


x(t)  ■  Ax(t)  +  Bu(t)  +  5*(t)  (1) 

where  u(t)  represents  the  control  Input  generated  by  the  control  law. 


Assuring  that  the  glldeslope  (and  of  coursa  localizer)  has  b««n 
captured,  the  landing  approach  (not  considering  flare,  touchdowi  and 
rollout)  bee  ones  a  tracking  problem  with  respect  to  the  deviation  froa 
the  reference  glides  lope  track  and  with  respect  to  the  nominal  or  trln 
values  of  the  other  aircraft  notion  variables.  For  a  fixed  reference 
glldeslope  track  angle,  this  Is  a  regulator  problea  free  a  control  point 
of  view.  Thus  the  eodel  for  the  flight  control  systaa  Is  developed  froa 
the  control  theor;  relating  to  the  optlaal  regulator  problea. 

As  a  brief  review  of  the  optlaal  control  regulator  problea,  con- 
aider  the  linear  systaa 

*(t)  •  A*(t)  ♦  bu(t)  ♦  C*(t)  (1) 

where  the  observations  are  of  the  fora 

y(t)  •  H«(t)  ♦  n(t)  (2) 

where  n  Is  a  zero  neen,  Gaussian  aaplltude,  white  noise  process.  Let 
J  he  a  quadratic  functional  In  state  and  control  deflnsd  hr 

J(u)  •  Lla  C  j *(T)  Qa(T)  ♦  u(T)Mhi(T)  |  (3) 

idiere  Q  and  A  are  syanetrlc,  Q  Is  non-negative  definite,  and  I  Is  posi¬ 
tive  deflnl*  .  It  Is  well  known  that  the  control  lew  defined  on  the 
observation  y,  which  alnlalze  the  functional  (3),  subject  to  the  dif¬ 
ferential  constraint  (1)»1s  of  the  fora  (Aef  5) 

u(t)  •  (4) 

where  F  Is  the  optlaal  feedback  for  the  deternlnlstlc  problea  (l.e., 
C*-0)  and  Is  given  by  the  equation 

F  •  -«-*rP  (I) 
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and  P  Is  the  solution  to  tho  nonlinear  matrix  Rlcattl  equation 

AP  ♦  PA'  ♦  Q  -  PBR^B'P  ■  0  (6) 

In  Cq  (4),  x*  is  the  "best  estimate"  of  x  and  Is  the  output  of  a  system 
of  differential  equations  (called  a  Kalman  filter)  **h1ch  have  as  their 
input  the  observations  y. 

Based  on  this  familiar  optimal  control  result  for  the  so-celled 
linear,  quadratic,  Gaussian  problem,  the  flight  control  system  Is  modeled 
by  the  cont^l  law 

u(t)  •  Fi(t)  (7) 

where  F  Is  the  optimal  feedback  matrix  for  a  quadratic  cost  function  of 
th«  form  given  by  £q  (3),  and  x(t)  <s  a  "sub-optimal*  estimate  of  the 
state  s,  and  is  described  In  Ref  1. 

The  quadratic  cost  funct'on  used  to  define  F  Is  of  the  form 

J»(u)  •  a\  ♦  k»o!  ♦  k,o.  (•) 

■  *C,  *c, 

mnere  o*  Is  the  rms  glide  slope  deviations  from  the  reference  gl  ides  lope 
treckjoj  ,  1  *1,2,  are  the  rms  control  Inputs.  For  the  DC -4  these 
control  Inputs  are  commanded  elevator  position  and  commanded  engine  r*n. 
In  the  case  where  direct  lift  control  Is  considered,  these  control  Inputs 
are  commanled  el  •rater  position  aiv<  commanded  direct  lift.  The  weighting 
coefficients,  A»  and  k«,  were  selected  to  mlnlmlie  the  prohablllty  of  e 
missed  eppreech  (PUR)  subject  to  prospect  fled  limits  on  the  rms  control 
ectlvlty.  le  this  manner,  the  structure  of  the  solution  to  the  llneer, 
quedretlc,  Ceusslen  optimal  control  problmi  uos  used  to  define  the  feed- 
beck  galea,  end  the  cost  functions!  (!)  wet  determined  so  es  to  opt  lot  to 
tho  lending  approach  performance  smile  ecceeetleg  for  the  control 


authorities  that  are  Imposed  on  the  flight  control  system.  The  value  of 
the  structure  of  the  solution  to  the  linear,  quadratic,  Gaussian  problem 
Is  that  Eqs  (5)  and  (6)  for  the  optimal  feedback  matrix  F  can  be  easily 
and  rapidly  solved  via  modern  digital  computer  techniques.  Thus  the 
problem  Is  reduced  to  solving  for  kt  and  kt  of  Eq  (3)  which  minimizes 
PHA  subject  to  specified  rms  constraints  on  the  control  activity.  This 
was  done  metrically  by  a  direct  search  minimization  technique  called 
“pattern  search.”  The  details  of  this  minimization  technique  and  the 
digital  computer  Implementation  are  given  In  Ref  4. 

As  previously  mentioned,  a  “sub-optimal*  state  estimate  x  Is 
used  In  the  control  law  formulation  given  by  Eq  (7).  The  form  of  the 
“sub-optimal*  state  estimator  and  the  reasons  for  using  the  “sub-optimal* 
stele  estimate  Instead  of  an  optimal  estimate  of  the  state  In  the  control 
definition  are  given  In  Refs  1  and  4.  This  aspect  of  the  optimal  control 
law  Is  not  particularly  pertlrent  to  this  study  tinea  It  Is  the  optimal 
feedback  gains  which  we  really  want.  The  Interested  reader  can  find 
complete  coverage  in  Refs  1  and  4. 

Since  the  control  lew  given  by  .r4  (7J  was  developed  for  the  case 
whore  disturbance  and  meiwreorni  noise  are  zero  mean  and  Gaussian,  u 
does  not  directly  deal  with  cases  where  steady  winds  and  wind  shears  ars 
considered.  This  approach  was  taken  to  avoid  the  Intractability  of 
dealing  with  a  probabilistic  description  of  these  types  of  disturbance. 

It  would  else  be  unreal Istlc  to  prescribe  o  flight  control  low  which  took 
into  account  a  given  deterministic  wind  directly  because  they  are  In  fact 
random  from  day  to  day  and  from  one  geographic  location  to  another.  It 
Is  possible,  however,  to  develop  a  control  law  which  suppresses  the  effect 
tf  steady  winds  and  wind  shear  &*>  the  deviation  from  glide  slepe.  This 


Is  done  by  Including  a  state  In  the  system  equations,  Eq  (1),  which  Is 
the  Integral  of  glldeslope  deviation.  In  this  way,  the  control  law  given 
by  Eq  (7)  Incorporates  an  Integral  control  which  suppresses  steady  errors 
In  the  glldeslope  deviation  due  to  steady  winds  and  compensates  for  the 
wind  shears.  This  approach  was  taken  for  the  DC-8  using  Integral  feedback 
In  glldeslope  tracking.  Recall,  however,  that  steady  winds  and  wind  shear 
are  taken  as  zero  In  this  stu4y. 

The  resulting  automatic  flight  control  system  model  used  In  the 
data  rate  analysis  Is  shown  schematically  In  Fig.  4.  The  discrete  Kalman 
filter  and  the  Intersample  extrapolator  shown  in  the  figure  are  used  to 
derive  the  state  estimate  x  and  are  described  In  Ref  1.  The  portion  of 
the  control  system  model  that  Is  of  particular  Interest  In  this  stu<1y  Is 
the  optimal  feedback  gains. 
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Fig.  4.  Block  Diagram  of  the  Landing  Approach  Model  with  the  Optimal  AFCS 


IV.  Standard 


In  this  chapter  a  classically  designed,  advanced  automatic  flight 
control  system  Is  described.  This  AFCS  Is  taken  from  Ref  2  and  It  Is 
used  as  a  standard  In  this  study.  This  standard  served  to  validate  the 
computational  techniques  of  the  digital  computer  program  and  was  used  to 
compute  baseline  values  for  such  figures  of  merit  as  the  probability  of 
missed  approach  (PMA),  rms  glldeslope  deviation  (04),  rms  pitch  angle 
deviation  (09),  rms  elevator  angle  deviation  (o$e)  and  the  rms  elevator 
actuator  rate  (a&e).  These  baseline  figures  are  used  as  a  measure  of 
"goodness"  with  which  to  compare  the  sub-optimal  AFCS  design.  Included 
In  this  chapter  are  the  feedback  control  law  for  the  standard,  the  state 
equation  formulation  for  the  standard,  and  a  comparison  of  numerical 
results  from  Refs  2,  6,  and  this  study  for  the  various  figures  of  merit. 

The  standard  automatic  flight  control  system  was  designed  using 
"classical"  multi-loop  control  techniques,  and  the  design  philosophy  as 
well  as  the  control  system  Is  described  In  Ref  2.  A  block  diagram  of 
the  standard  system  Is  shown  In  Fig.  5.  The  standard  does  not  represent 
an  existing  system,  but  rather  an  advanced,  high  performance  automatic 
pilot  and  approach  coupler  of  the  type  that  would  be  required  for  suc¬ 
cessful  Category  II  landing  approach  operations  In  a  moderate-to-severe 
turbulence  environment. 

The  Inner  loop  of  Fig.  5  Is  used  to  feed  back  pitch  and  pitch  rate. 
Pitch  rate  6  Is  fed  back  for  short-period  damping  and  to  extend  the  path¬ 
following  bandwidth.  A  conventional  feedback  of  pitch  attitude  6  to 
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Block  Diagram  of  Automatic  Glidepath  Displacement  Control  System  of  the  Standard 


achieve  short-period  stiffness  together  with  path  damping  Is  modified 
by  a  washout  with  a  comparatively  small  time  constant  T^.  This  has 
the  effect  of  retaining  the  short-period  attitude  stiffness,  but  It 
trades  path  damping  for  a  reduced  glide  path  displacement  due  to  normal 
gusts.  The  gains  and  time  constants  used  In  the  Inner  loop  closure  are 
given  In  Table  III. 

TABLE  III 

Selected  Gains  and  Time  Constants  for  the  AFCS  of  Ref  2 

Pitch  Rate  and  Attitude  Stability  Augmentation 
1/Ta  ■  15.0  rad/sec 
1/Two  “  0,7  ra£,/sec 

■  -2.0  sec 

Kq  -  -2.0 

Path-Following  Regulation  and  Control 
1/Tf  -  2.0  rad/sec 

■  -0.00867  rad/ft 
K j  -  -0.0013  1/ft-sec 
Kj  ■  -0.0256  sec/ft 


The  outer  loop  provides  the  feedback  and  filtering  of  glldeslope 
displacement  Information  and  Is  referred  to  as  a  path -foil owing  loop. 
Since  the  pitch  attitude  6  feedback  Is  modified  by  a  washout,  the  air¬ 
craft  has  no  absolute  attitude  reference.  This  dictates  the  need  for  a 
high  quality  path  damping  signal.  In  practice,  this  could  be  the  derived 
beam  rate  d,  the  Incremental  altitude  rate  h,  or  the  result  of  a  comple¬ 
mentary  filtering  which  could  Include  output  of  a  normal  accelerometer. 
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Tht  analysis  presentee  In  Ref  2  assumed  thet  In  e  high-performance  system, 
end  on  •  reUtlvely  shallow  glldepath,  such  es  -3.0  deg,  the  complementary 
filtering  c«n  be  performed  In  such  e  way  thet,  over  the  bandwidth  of  the 
peth-fol lowing  loop,  e  pure  or  noise-free  d  signal  can  be  provided. 

A  been  displacement  signal  Is  required  for  path  acquisition  and 
stiffness.  To  this  Is  added  the  Integral  of  the  bean  displacement  to 
keep  the  aircraft  on  the  reference  glldepath  In  the  presence  of  a  head¬ 
wind  or  a  long  wavelength  updraft.  The  gain  on  the  Integral  term  Is 
limited  by  considerations  of  path-following  stability,  so  that  Its 
effectiveness  Is  only  felt  In  regulating  against,  at  most,  slowly  changing 
winds.  Both  the  beam  displacement  and  Integral  of  beam  dlsplacwent  are 
shown.  In  Fig.  5,  to  be  filtered  by  a  low-pass  filter  with  a  time  constant 
Tf.  This  filter  Is  representative  of  the  combined  Impedance  of  the  fil¬ 
ter  capacitor  and  the  receiver  conventionally  used  In  the  VHF-UHF  ILS. 
Alternatively,  Tf  can  be  taken  to  closely  approximate  the  combined 
characteristics  of  a  receiver  boxcar  hold  and  ripple  filter  such  as 
might  be  used  In  connection  with  a  microwave  scanning  beam  system.  The 
numerical  values  for  the  gains  and  time  constants  used  In  the  outer  loop 
closure  are  given  In  Table  III. 

The  feedback  control  law  for  the  system  of  Fig.  5  Is  arrived  ct  by 
closing  both  feedback  loops  and  assuming  dc»0.  The  control  law  6Cc  Is 
given  by 

6ec  ■  -  Y g®  *  0C 
The  first  term  Yg  Is  given  by 


Y 


6 


S  +  l/T*o 


♦ 


K,. 


K$»(s  ♦  1/Te) 
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where 


-L.fi.-L 

te  *i  T« 


Tht  second  torn  ec  Is  given  by 


where 


®c  *  -  V 


*d 


*d  **d» 
*(S  ♦  1/Tf) 


and 


Ka(s  ♦  1/Td  )(s  ♦  1/Td  )(s  ♦  1/Td  ) 
S{*  ♦  1/Tf) 


J_ 

Tf 


Expressing  tht  control  law  aquations  In  state-variable  for*  requires 
a  little  "fiddling  around"  and  the  establishment  of  at  least  three  dusryr 
states.  For  the  benefit  of  the  Inquisitive  reader  the  equations  used 

are 

ua$  ■  xuuas  ♦  V  -  ge  cos  y.  ♦  «ugug  -  X*wg  ♦  X6<,de  ♦ 
w  ■  ZuUi$  ♦  Z*w  ♦  U.q  -  ge  sin  yi  -  Itfq  ♦  Z«e«e  ♦  Zs^dth 
q  -  M„uts  ♦  M*w  ♦  M<,q  -  M*wg  ♦  M6#d#  ♦  M^d^, 

9  ■  q 

d  ■  U|6  -  w 
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4#  .  -  J-  4#  ♦  -L  4^ 

U  Tl 

if  •  •  -J-  $f  ♦  Kgq 


if  m  -  -L  df  -  6  -  y, 


#•  •  *f 


where  6f,  df ,  And  y,  if*  dua*y  itotit  representing  the  output  of  the 
washout,  the  output  of  tht  glldaslope  receiver  filter.  And  tho  Integral 
of  tho  fUterod  glldoslopo  deviation.  Tho  steto  y»  represents  tho  Addi¬ 
tive  nolso  pieced  on  tho  glldoslopo  dovlotlon  d.  Tho  gains,  tint  con¬ 
stants,  and  stability  derivatives  art  given  In  Tables  II  and  III. 


Comparison  of  Results 

Tho  nunerlcal  results  of  Ref  2  Mere  validated  by  tw>  methods:  one, 
using  tho  digital  computer  program  that  mis  used  In  this  study  to  compute 
tho  numerical  results,  and  tM>,  an  analog  simulation  technique  reported 
In  Ref  6.  The  comparison  Is  shoun  In  Table  IV. 

The  figures  In  Table  IV  shoe  that  exact  agreement  mas  not  achieved 
between  the  data  of  Ref  2  and  the  digital  solution.  Also  note  the  results 
of  the  analog  simulation  of  Ref  6  do  not  agree  exactly  with  those  of  Ref  2. 
The  digital  computer  results  and  those  of  Ref  6  tend  to  vary  In  the  same 
direction  when  compared  to  the  results  of  Ref  2.  The  root  sum  squared 
disturbance  correlated  error  figures  show  fairly  good  agreement  between 
all  three  studies,  and  the  values  for  04  and  09  of  Ref  2  are  used  as  a 
basis  for  comparison  In  this  study.  A  value  for  the  ms  elevator  rate 
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TABLE  nr 

Standard  Deviations  of  Systea  Response  to  Stochastic  Disturbances 
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1«  not  glvfii  In  Mf  2;  thtrtfort.  •  voluo  of  .1692  r«d/»oc,  «s  coaputod 
0/  tho  digital  progroa,  will  bo  usod  for  coopering  control  octlvltlot. 
Tho  rotponto  of  tho  'lUndird*  tyttoa  to  •  pitch  roto  lapulso  It  thorn 
la  Fig.  12  of  ChtpUr  VII. 
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V.  Developing  the  Practical  System 


The  aeln  objective  of  this  study  Is  to  determine  If  «  practical 
control  system  could  be  developed  froei  •  full  optimal  system  through  soee 
logical  procedure.  The  step  by  step  procedure  Is  presented  for  two  casts 
using  the  equations  of  notion  for  the  DC-8  aircraft  as  given  In  Chapter  II. 
In  the  first  case,  airspeed  was  assumed  to  be  Matured  and  available  for 
the  control  law,  while  the  second  case  did  not  use  airspeed.  Before 
starting,  a  brief  review  of  the  complete  procedure  Is  given. 

As  explained  In  Chapter  II,  the  optlMl  control  law  Is  developed 
to  mlnlmlxe  a  given  cost  functional  subject  to  certain  constraints.  In 
order  to  arrive  at  this  control  law,  It  Is  necessary  to  specify  the  Mature- 
Mnts  which  are  assuMd  to  be  available.  As  previously  Mntloned,  tone  of 
these  MasureMnts  are  associated  with  the  flight  control  system  and  are 
available  continuously,  free  of  noise.  Those  MasureMnts  considered  here 

Include  sensed  pitch  attitude,  pitch  rate,  and  noraal  acceleration.  These 

• 

are  coanon  to  both  cases  presented  In  this  chapter.  The  regaining  Masure¬ 
Mnts  are  the  glldeslope  deviation,  the  associated  errors  as  Mntloned  In 
Chapter  II,  and  for  the  first  case,  airspeed  with  noise  and  MasureMnt 
errors.  The  Integral  of  glldeslope  deviation  Is  derived  from  the  glide- 
slope  Masurement. 

Once  given  the  desired  measurements,  the  optimal  control  law  can 
be  solved  to  obtain  the  optimal  feedback  gains.  The  feedback  gains  Include 
a  set  of  gains  for  an  autothrottle.  Since  the  approach  flight  condition 
Is  at  a  trim  speed  above  the  speed  for  trim  response  reversal,  both  glide- 
path  tracking  and  speed  regulation  can  be  achieved  with  elevator  control 
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•loot  (Rtf  1).  Si  net  autothrottle  ft  not  needed,  Uit  first  step  la  to 
remove  those  feedbacks  (l.t.,  rtro  tht  appropriate  gains).  The  procedure 
continues,  removing  feedbacks,  one  at  a  tint,  until  only  those  desired 
feedbacks  in  ltft.  Of  those  flntl  feedbacks,  glfdeslope  deviation  and 
•Irspttd  nutt  have  measurement  noise  added  «nd  then  bt  filtered.  Tht 
filtered  measurements  trt  used  In  tht  flntl  control  It*  fine#  they  ftp re¬ 
sent  tht  Infornttlon  it  It  would  really  bt  available  for  use  In  tht  AfCS. 

If  tht  fttdbtck  control  low  It  simplified  further  by  removing  feed* 
becks  on  measurements  assumed  tvt titbit  In  tht  opttmiiatton  procett,  tht 
resulting  control  low  will  bt  tlthtr  unstable,  seriously  degrade  approach 
performance,  result  In  high  tc  tut  tor  rates,  or  result  In  poor  pitch  Attitude 
chtrtcUrlsttcs.  Tht  results  of  deleting  tht  presptclfltd  measurements  ere 
given  In  Table  V  end  reinforce  tht  origlntl  Hypothesis  (l.t.,  the  nets ure- 
ments  stifled  for  optlmitation  nwst  bt  retained  by  tht  practical  systml. 

All  dtU  presented  In  this  ptptr  tssuned  •  sctiwilng  been  system 
•s  described  In  Rtf  1.  A  stapled  dttt  re tt  of  t  senpUs/tec  wes  used  for 
til  casts.  This  figure  wes  found  to  bt  on  acceptable  dtu  rate  for  an 
AFCS  Incorporating  normal  acceleration  Information  in  tht  control  law. 

Control  law  with  Airspeed  Feedback 

In  the  first  cast,  airspeed  information  Is  assumed  to  bt  available 
In  formulating  the  control  law.  Tht  procedure  Is  described  a  step  at  a 
time  and  the  results  are  ecplalned  at  each  step.  Tht  numerical  data  for 
each  sttp  are  given  In  Tablt  V. 

Step  1.  Tht  equations  are  solved  using  the  full  state  feedbacks 

from  tht  optimisation  routine  (data  rate  analysis  program). 
This  gives  the  performance  attainable  with  full  Uata 
feedback. 
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Sup  3.  The  opt  imitation  procure  two  control  mt 

•levator  control  law  and  a  throttle  control  tow.  AuU- 
.  throttle  wot  not  used  fn  Ref  l  for  the  Standard  and  l| 

not  M»ded,  to  UO  «wlO thro tile  fetdMMS  or*  deleted. 

Stop  ).  (levator  pot  it  <  on  ft  fed  Met  in  the  fwtt  optimal  system 
The  standard  does  not  uta  tltvaUr  position  to  this  feed- 
Met  it  dtltted. 

Sup  4.  Tht  optima)  feedMek  gains  inctwdt  feedMek  of  Uo  wind 
gust  intensities,  Mind  gusts  cannot  M  measured  easily 
SO  feeding  them  Met  If  not  practical.  Thesa  Mint  or* 
tot  U  roro. 

At  this  point,  in#  control  tow  consists  of  only  mo  dotirod  food* 
Mens;  oirtpood  w#|,  normal  voloclty  w,  pitch  attitude  *.  pitch  rou  Q, 
glides  lope  deviation  d,  glides lope  deviation  rtu  4,  ond  mo  integral  of 
glldetlope  dev lotion  Jif.  To  verify  tho  constraint  on  tha  typothesit  of 
CMpur  I.  tho  protpocified  foodMett  oro  now  individually  zeroed. 

Stop  S.  Zero  tho  feedMek  on  oirtpood  u^,  the  PKA  increases. 

Stop  €.  Zero  the  foodMek  on  no  nut  velocity  w.  ond  tho  fttA 
lncfT«tet. 

Stop  7.  Zero  the  foodMek  on  ptuh  attitude  0,  ond  tho  system 
It  wntublo. 

Stop  6.  Zero  the  foodMek  on  pitch  rou  q.  This  stop  desorvos 
spec  it  t  cowont.  The  previous  stops  showod  a  definite 
dotrioonul  effect  on  tho  porforwnce  while  mis  step  Ms 
improved  the  approach  performance  PKk.  hoticc,  however, 
tMt  tht  ms  pitch  ottitudt  doviotion  Ms  increased  os  Ms 
the  ms  ectuetcr  rote,  Removing  pitch  rate  feedback 
decreases  the  response  tine  but  causes  the  system  to  Mve 
poor  short-period  damping.  The  light  cooping  is  more 
evident  when  the  system  response  to  a  pitch  rou  impulse 
is  plotted  es  in  Fig.  6.  Inasmuch  at  MS  is  computed 
usino  the  ms  gif  deslope  deviation,  these  high-frequency 
oscillations,  Mvfng  •  low  ms  cooponent,  contribute  littlw 
to  the  ms  alideslope  deviation  and  FNA.  Iftilla  tht  PNA 
is  accepubla,  the  hi  ah -frequency  oscillations  produce 
uMccepUble  rlda  qualities  and  unaccepUbla  pitch  attitude 
characteristics. 

Step  9.  Zero  tha  feedteck  on  glides  lope  deviation  d,  and  tht 
system  is  unsUblt. 
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o  oo  ej.oo  40.oo  60.00 


TIME  (SEC.  ) 


f\%.  6.  Sjnu*  iHpOAM  tfthowt  Mun  llU  fllOiel 


Stop  10.  Zero  (ty  Oft  tty  integrel  of  glidttlope  depletion  /d 

one  tty  control  tyttoo  will  no  longer  cooper  tote  for  J 
it 0*0/  viMi. 

At  type  lAet  I /fd  to  Ctypter  1,  Uw  roopfil  of  Of/  p  ret peel  fled  feed- 
boa  rote)  It  to  oo  t*te  tuple  control  tyttoo.  Aw  returning  to  tty  control 
tyttoo  oo  It  mi  configured  to  Stop  4,  tty  glidttlope  demotion  end  elr- 
tpeed  foo  Ami  oust  bo  oodt  u  rtsooblo  tty  *rto»  world*  ooofurooonts. 

Stop  11.  Tty  gltdfflopo  deration  oootwrcMirftt  oolft  y,  ond  orrort 
twe*  Of  fipUlnod  to  Diopter  It  or*  conftdorod  first.  Tty 
ootfo  it  oddod  to  tty  glldetlop*  depletion  ond  o  now  tuto 
tf  introduced  to  roprtfoot  tty  filleted  el  Ides lope  deration 
oooteteteftt  df.  Tty  filter  if  ftyun  to  rig.  7.  Tty  tint 
comUM  If  not  cbosoo  to  ototofto  tty  fWA.  Tty  reluo 
selected  Ml  0.1  toe.  It  not  Of  flood  tty  df  tt  olto  tff«d 
to  dor  let  tty  totofrol  of  d  foodboa. 


t _ .•/ 

Pi 

• 

S _ - 

1 

•f^ 

- - - 

Iff  •  1 

ftg.  7.  Oltdfflopo  Deration  Wrote rtoont  Model 


Stop  17.  Including  eirtpred  In  tty  control  low  requires  ttyt  tty 

problem  oifoctotfd  writli  oootering  eirtpred  bo  contidtred. 
Ttyro  tt  no  I  it  on  tty  o—feroroM.  but  ooro  ioportont,, 
ttyrt  (tint  significant  tint  log  In  IM  HMtercwtnt 
cowvod  by  tty  ooctyntfo  used  to  mil  tbr  oootwtenont. 

Airs  prod  OMteroomt  dewlcet  currently  tn  eto  incur  o 
tioo  drloy  In  tty  order  of  O.S  toe.  It  it  rtotonoblt  to 
ottwor  o  log  to  tontod  oirtptod  of  O.S  toe.*  Tty  outturn- 
ooot  drvico  wot  opprootootod  by  o  first -order  log.  It  wot 


•Soffit tod  by  Aooold  Andrrsoo,  Air  forco  flight  Oynooict  Uborttory, 
Mrtgbl-Petlerton  All,  Obto. 


TABLE  V 


Numerical  Results  of  Sub-optimal  Control  System  with  Airspeed  Feedback 


Step 

PMA 

3o  Value  of 

Elevator  Actuator  Rate 
(rad/sec) 

lo  Value  of 
Pitch  Angle 
(rad) 

1 

.02084 

.16572 

.02119 

2 

.02090 

.16572 

.02153 

3 

.01719 

.16974 

.02178 

4 

.01074 

.23367 

.02436 

5 

.19926 

.04219 

.02097 

6 

.24631 

.24463 

.01781 

7 

unstable 

8 

.00187 

1 

.27216 

.03177 

9 

unstable 

10 

.00468 

.23137 

.02374 

11 

.01235 

.23932 

.02422 

12 

.03980 

.10700 

.02518 

Standard 

.0366 

.168 

.0287 
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assumed  that  most  of  the  sensor  noise  present  would  be 
filtered  out  by  the  device  Itself.  The  measurement 
model  Is  shown  In  Fig.  8,  where  Ti-0.5  sec.  The 
measured  airspeed  is  called  um.  The  system  performance 
Is  extremely  sensitive  to  the  airspeed  measurement.  A 
delay  of  1.0  sec  was  found  unacceptable  as  the  system 
performance  deteriorated  greatly  (PMA». 07304). 


uas 

1 

um 

Tls  +  1 

Fig.  8.  Airspeed  Measurement  Model 


The  results  of  Step  12  represent  the  final  system.  The  numerical 
results  are  compared  with  those  for  the  standard  in  Chapter  VII.  The 
system  response  to  a  pitch  rate  Impulse  of  1  deg/sec  Is  shown  In  Chapter 
VII,  Fig.  13. 

Control  Law  Without  Airspeed  Feedback 

This  case  Is  developed  much  as  the  previous  case  except  airspeed 
Is  not  considered  as  one  of  the  measurements  available  In  the  development 
of  the  optimal  feedback  control  law.  Developing  the  control  law  with  the 
same  measurements  as  those  used  In  Ref  2  permits  the  comparison  of  a  clas¬ 
sically  designed  advanced  system  with  the  realizable  sub-optimal  system 
where  the  same  feedbacks  are  assumed  to  be  available.  The  numerical 
results  for  this  case  are  found  In  Table  VI.  The  step-by-step  procedure 
Is  as  follows: 

Step  1.  As  In  the  first  case,  the  Initial  step  Is  to  use  the 
full  state  feedback. 


41 


Step  i.  delete  all  feedbacks  to  the  throttle. 

Step  3.  Zero  the  elevator  position  feedback. 

Step  4.  Zero  the  wind  gust  Intensity  feedbacks. 

Step  5.  Zero  the  airspeed  feedback.  It  should  be  pointed  out 
here  that  even  though  airspeed  was  not  a  prespecified 
feedback,  airspeed  feedback  does  have  a  significant  effect 
on  approach  performance  and  control  activity. 

The  control  law  now  consists  of  just  the  desired  (prespecified) 
feedback  gains,  and  measurement  noise  or  filtering  Is  not  yet  considered. 

As  In  the  previous  case,  these  desired  feedback  gains  are  now  Individually 
set  to  zero,  one  at  a  time,  to  validate  the  constraint  on  the  hypothesis 
In  Chapter  I  and  to  show  the  resultant  effect. 

Step  6.  Zero  the  normal  velocity  feedback. 

Step  7.  Zero  the  pitch  attitude  feedback. 

Step  8.  Zero  the  pitch  rate  feedback. 

Step  9.  Zero  glldeslope  deviation  feedback. 

Step  10.  Zero  the  Integral  of  glldeslope  deviation  feedback.  The 
improvement  In  this  step  Is  expected  since  Integral  feed- 
jacks  have  a  destabilizing  effect.  The  cost  of  this 
Improvement  Is  the  loss  of  compensation  for  steady  winds. 

With  fewer  measurements  available  for  the  control  law,  the  gains 
are  higher  and  the  result  of  removing  one  of  these  feedbacks  Is  much  more 
pronounced.  This  Is  apparent  when  comparing  the  effects  of  Steps  6  through 
10  of  Table  VI  with  the  corresponding  steps  of  Table  V.  Now  returning  to 
the  control  system  as  It  was  configured  In  Step  5,  the  glldeslope  deviation 
must  be  made  to  resemble  the  actual  measurements. 

Step  11.  The  glldeslope  deviation  Is  corrupted  with  noise  y2,  and 
filtered  as  In  the  previous  case.  Originally  the  filter 
In  Fig.  7  was  used;  however,  the  effect  of  filtering 
glldeslope  deviation  was  a  reduced  elevator  actuator  rate 
and  consequently  a  poor  PMA.  The  filter  of  Fig.  7  was 
modified  by  Including  a  gain  6  in  the  numerator  to 
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compensate  for  the  effect  of  filtering.  This  gain  was 
varied  to  achieve  the  minimum  PMA.  This  filter  is  shown 
in  Fig.  9,  where  Tf  is  0.5  sec  and  G  is  1.27. 


ya 

* .  A.  r 

G 

df 

Tfs  +  1 

Fig.  9.  Glideslope  Deviation  Measurement  Model  with  Increased  Gain 


Step  12.  The  control  system  of  Step  11  has  a  much  poorer  per¬ 
formance  than  the  standard.  The  standard  achieved  a 
PMA  of  .0366  while  the  sub-optimal  control  system,  at 
this  point,  results  In  a  PMA  of  .07316.  A  closer  look 
at  the  data  reveals  that  one  cause  for  the  poorer  per¬ 
formance  Is  a  much  lower  rms  elevator  actuator  rate: 
.08058  rad/sec  for  the  sub-optimal  as  opposed  to  .16779 
rad/sec  for  the  standard.  The  logical  "fix"  is  to 
Increase  the  rms  elevator  actuator  rate  in  the  sub- 
optimal  system.  To  get  the  higher  actuator  rate,  the 
optimization  routine  was  rerun  with  relaxed  constraints 
on  actuator  activity.  This  resulted  in  higher  feedback 
gains  and  In  the  final  system,  higher  rms  actuator  rates. 
It  Is  Important  to  point  out  the  ease  and  quickness  with 
which  this  Improvement  was  achieved.  Because  the  equa¬ 
tions  are  mechanized  for  computer  solution,  only  one 
card  needed  to  be  changed  and  the  entire  revision  was 
completed  In  less  than  a  day.  The  final  system  response 
to  a  pitch  rate  impulse  is  shown  in  Chapter  VII,  Fig.  14. 

The  final  sub-optimal  control  law  results  In  an  oscillation  when 
the  system  1r  perturbed  by  a  pitch  rate  impulse  as  shown  by  Fig.  13  of 
Chapter  VII.  A  classically  designed  system  (Ref  7)  similar  to  the 
standard  but  without  the  washout  also  experienced  those  oscillations. 

If  these  oscillations  are  considered  undesirable,  the  pitch  attitude 
feedback  can  be  modified  by  a  rashout.  This  has  the  effect  of  retaining 
the  short-period  attitude  stiffness,  but  It  trades  path  damping  for  a 
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reduced  response  or  glide  path  displacement  due  to  normal  gusts  or  pitch 
rate  impulses. 

The  realizable  systems  described  in  this  chapter  will  be  compared 
to  the  advanced,  classically  designed  standard  in  Chapter  VII.  The 
numerical  values  for  the  feedback  gains  are  given  In  Appendix  A.  The 
directness  and  simplicity  of  this  technique  should  be  quite  evident  at 
this  point. 


44 


TABLE  VI 

Numerical  Results  of  Sub-optimal  Control  System  Without  Airspeed  Feedback 

(6  samples/sec) 


Step 

PMA 

3o  Value  of 

Elevator  Actuator  Rate 
(rad/sec) 

lo  Value  of 
Pitch  Angle 
(rad) 

1 

.00131 

.26750 

.02443 

2 

.00131 

.26750 

.02443 

3 

.00100 

.27791 

.02469 

4 

.00136 

.34038 

.02633 

5 

.05172 

.07137 

.02250 

6 

.31705 

.14866 

.01940 

7 

unstable 

8 

unstable 

9 

.79171 

.06014 

1 

.03216 

10 

.05173 

.12369 

.01895 

11 

.07316 

.08058 

.02510 

12 

.02154 

.15095 

.02680 

Standard 

.0366 

.168 

.0287 
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VI.  Application  to  Direct  Lift  Control 


In  this  chapter,  the  technique  described  in  Chapter  V  is  used  to 
develop  realizable  control  laws  for  a  hypothetical  aircraft  which  has 
the  capability  of  using  direct  lift  control.  Direct  lift  control  has 
attracted  considerable  attention,  particularly  for  possible  use  on  STOL 
aircraft.  The  primary  purpose  for  considering  DLC  is  to  conceptually 
show  the  significant  improvement  which  can  be  achieved  in  glideslope 
tracking  and  to  demonstrate  the  versatility  of  this  design  technique. 
Direct  lift  control  represents  a  multiple  input,  multiple  output  system 
which,  for  classical  control  techniques,  is  at  best  extremely  difficult. 
The  optimal  approach  handles  this  case  with  comparative  ease. 

This  chapter  describes  how  DLC  was  conceptually  considered  for  the* 
DC-8.  The  modified  equations  of  motion  and  actuator  equations  are  given. 
Two  cases  are  presented  in  detail:  the  first  case  includes  airspeed 
measurements  in  the  control  law,  while  the  second  case  does  not.  The 
procedure  follows  closely  that  described  In  the  previous  chapter.  The 
numerical  results  for  each  step  of  the  procedure  are  in  Table  VIII.  The 
system  response  to  a  pitch  rate  Impulse  Is  found  In  Chapter  VII,  Figs.  15 
and  16.  The  numerical  values  of  the  optimal  feedback  gains  are  presented 
In  Appendix  A. 

Aircraft  Model  with  Direct  Lift  Control 

The  Investigation  of  DLC  was  performed  on  a  hypothetical  DC-8 
incorporating  Idealized  direct  lift  control.  The  model  Is  taken  from 
Ref  1,  where  It  was  assumed  that  Ideal  DLC  would  generate  normal  forces 
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on  the  aircraft  without  significant  changes  in  pitch  attitude.  The  basic 
aircraft  equations  were  adapted  to  include  direct  lift  control  by  the 
introduction  of  a  new  control  input  vector 


u  = 


°ec 


where  6<jc  is  the  commanded  input  to  the  direct  lift  control  surface. 
Thus  the  throttle  control  was  replaced  by  a  hypothetical  direct  lift 
control . 

The  equations  of  motion  in  Chapter  II  were  modified  by  replacing 
the  thrust  stability  derivatives  (X,$  (  . ,  Z^,  M^)  with  a  hypothetical 
set  of  direct  lift  stability  derivatives  (X5  ,  Z^,  M^).  The  equations 
for  the  hypothetical  "direct  lift"  DC-8  are 

“as  E  xuuas  +  xww  -  go  cos  70  +  mUgug  -  Xwwg  +  X6e6e  +  Xfi^ 


W  «  ZuUas  *  Zww  -  ge  sin  y0  «■  boq  -  !wwg  +  +  2$^ 


q  -  Muuas  +  Mww  +  Mqq  -  Mwwg  +  M«e6c  +  M(Sd6d 
0  *  q 

d  *  Uo6  -  w 

The  values  for  the  direct  lift  stability  derivatives  are  given  In  Table 
VII.  The  remaining  stability  derivatives  are  left  the  same  and  are  found 
In  Table  II.  The  disturbances  are  the  same  as  described  in  Chapter  II. 
The  differential  equations  for  the  elevator  actuator  and  direct  lift 
control  surface  are 
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60  ~  -  —  6p  ^  —  6 p . 

a  t  H  t  ec 
'e  'e 


6d  =  -  ~  6d  +  --  6d 
Td  Td  c 

where  6d  is  the  direct  lift  control  surface  deflection,  Td  and  Te  are  the 
control  lags  on  the  direct  lift  surface  and  the  elevator;  Td  and  Te  were 
both  taken  to  be  .06666  sec. 

TABLE  VII 

Direct  Lift  Stability  Derivatives  (Ref  1) 

X6d  =  0.0  ft/rad-sec2 

Z6d  =  -50.0  ft/rad-sec2 
Mfid  -  0.0  1 /rad-sec2 


The  numerical  value  for  76d  corresponds  to  a  control  effectiveness 
of  roughly  0.1  g  nonnal  acceleration  for  a  4°  surface  deflection.  The 
choice  of  X6d*M{d*0  causes  the  direct  lift  control  to  be  relatively 
uncoupled  from  the  G,  q,  and  uflii  equations.  This  Is  especially  true 
since  Mw  is  zero. 

Control  Law  with  Airspeed 

The  first  case  considered  for  the  hypothetical  aircraft  with  DLC 
Included  airspeed  feedback  in  the  control  law.  The  other  measurements 
used  In  the  optimization  process  were  normal  velocity,  pitch  attitude, 
pitch  rate,  and  sampled  g’ideslope  deviation.  The  procedure  follows 
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closely  that  used  in  the  previous  chapter.  There  are  two  control  laws 
of  interest  when  considering  DLC:  the  elevator  control  law  and  the  direct 
lift  control  law.  Therefore,  each  variable  has  two  feedback  gains.  The 
numerical  results  are  given  in  Table  VIII.  The  step-by-step  procedure 
is  as  follows: 

Step  1.  The  equations  are  solved  using  full  state  feedback  to 
evaluate  the  system  capability  using  the  full  set  of 
optimal  feedbacks. 

Step  2.  The  wind  gust  intensity  feedbacks  are  deleted  because 
wind  gusts  are  not  measured  in  a  realizable  flight 
control  system. 

Step  3.  Feedback  gains  on  the  position  of  both  the  elevator  and 
the  direct  lift  control  surface  are  set  to  zero.  This 
case  represents  the  prespecified  feedbacks  without  con¬ 
sidering  measurement  noise  or  filtering  of  the  measured 
signals . 

The  following  series  of  steps  were  used  to  verify  the  original 
hypothesis  of  Chapter  I  which  stated  that  removing  a  measurement  from  the 
control  law  which  was  assumed  present  for  the  original  optimization  results 
in  a  poor  control  system.  In  Table  VIII,  case  (a)  represents  the  result 
when  the  particular  feedback  was  deleted  from  the  elevator  control  law, 
while  case  (b)  is  the  result  when  the  some  feedback  was  deleted  from  the 
direct  lift  control  law. 

Step  4.  Zero  the  airspeed  feedback  gain. 

Step  5.  Zero  tie  normal  velocity  feedback  gain. 

Step  6.  Zero  the  pitch  attitude  feedback  gain. 

Step  7.  Zero  the  pitch  rate  feedback  gain.  Deleting  pitch  rate 

feedback  results  in  lower  response  times  but  higher  settling 
times.  Just  as  in  similar  cases  in  Chapter  V,  the  PMA  is 
satisfactory  because  the  high  frequency  oscillations  con¬ 
tribute  little  to  the  rms  glideslope  deviation.  The 
oscillations  cause  a  marked  increase  in  the  rms  pitch 
attitude.  Removing  pitch  rate  feedback  results  in  poor 
rice  qua  1 i ty. 
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-ABLE  VIII 

Numerical  Results  of  SuO-optlmal  Direct  Lift  Control  System 
Wltn  Airspeed  feedback  {Data  rate,  6  sarr.ples/sec) 


Step 

PKA 

3o  Value  of 
Elevator  Actuator  Rate 

(rad/sec) 

3o  Value  of 
Direct  Lift 
Actuator  Rate 
(rad/sec) 

lo  Value  of 
Pitch  Angle 

1 

.02428 

|  . 029C0 

.19461 

.01163 

2 

.00050 

.04361 

.27452 

.01187 

3 

.00035 

.04423 

.27898 

.01199 

4a 

.00596 

.00254 

.27924 

.01114 

b 

.16504 

.04454 

.03987 

.01250 

5a 

.00162 

.04546 

.27907 

.00997 

b 

.08778 

.04443 

.29341 

.01163 

6a 

.00146 

.04554 

.27968 

.01561 

b 

.12760 

.04463 

.29146 

.01615 

7a 

.00028 

.04622 

.27889 

.01268 

b 

.00000 

.04423 

.29391 

.01247 

8a 

.00053 

.04425 

.27889 

.01149 

b 

.08345 

.04420 

.27892 

.01181 

9 

.00047 

.04420 

.27862 

.01189 

10 

.00280 

.01709 

.11152 

.01208 

Standard 

.0366 

.168 

n/a 

.0287 

NOTE:  Case  a  represents  the  result  when  the  particular  feedback  was 
deleted  from  the  elevator  control  law,  while  Case  b  is  the 
result  when  the  same  feedback  was  deleted  from  the  direct  lift 
control  law. 
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Step  6.  Zero  u«?  5 1  ici*s 3ypc  deviation  fecdWt*  gain. 

The  results  of  trie  uu  f •  vf  steps  vallate  the  origin*)  hypothesis 
of  Chapter  1.  Now  going  u*» >  to  the  control  system  in  Step  J,  the  sampled 
gli deslope  deviation  and  airspeed  k!1 *  he  handled  as  in  the  first  cose  in 
the  previous  chapter. 

Step  9.  The  gliotslope  auction  ^as  assuscc  to  be  sinpled  at 
a  rate  of  C  samples /sec  and  to  nave  the  saa  errors  as 
described  in  Chapi cr  1 1  (i.e.,  broadband  additive  noise). 
The  noisy  .  oasure.  ent  infor.%»tion  has  fed  into  a  filter 
as  shewn  in  Fi«j.  7.  The  filter  tire  constant  was  pitted 
to  mini  i?c  th*.  probability  of  rissed  approach.  In  this 
case,  the  value  selected  >*as  '<■•0.2  sec.  The  filtered 
glides  lope  deviation  was  then  fed  back.  The  filtered 
glides  lope*  deviation  was  also  used  to  derive  integral  of 
glideslopc  dovi lu... 

Step  10.  The  airspeed  was  nandled  :n  the  same  manner  as  in  the 

first  case  in  C  aptcr  V  As  Iwiforc,  it  was  assumed  that 
any  noise  on  the  measured  airspeed  would  effectively  be 
filtered  by  me  mechanism  of  the  measurement  device.  It 
was  also  assured  that  the  reu  jrement  would  incur  a  time 
lag  due  to  the  mechanical  nature  of  the  measurement 
device.  Trie  measurement  model  is  shown  in  Fig.  8  of 
Chapter  V.  The  tin*  lag  T^  was  again  taken  to  be  0.5  sec. 

The  realizable  control  system  gives  a  very  low  PXA  (.00280).  This 
demonstrates  the  benefit  of  using  DLC.  The  rnss  direct  lift  control  sur¬ 
face  rate  dropped  more  then  half  when  the  lag  was  placed  on  the  airspeed 
measurement.  This  points  out  the  strong  reliance  on  good  airspeed  measure¬ 
ment  for  any  AFCS  which  considers  airspeed  as  a  prespecified  measurement 
In  developing  the  optimal  feedbacks. 


Control  Law  Without  Airspeed 

Prespecifying  airspeed  in  the  control  law  creates  several  problems. 

As  the  previous  case  demonstrates,  the  control  system  with  airspeed  feed¬ 
back  tends  to  be  very  dependent  on  this  feedback.  This  dependency  results 
In  the  need  for  nigr.  quality  sensors  which  will  create  a  minimum  time  lag 
in  the  measurement  and  also  filter  out  most  of  the  noise.  These  disadvantages 


fro*  10#  f*»r  dctcloplns  4ft  OHH41  Control  la*  without  a'rtpeed 

*t^W  l  1*4  procrAw**  ,.»f*iic'v  th»t  df  It*  second  cate  In  Chapter  V. 

?R#  numerical  results  «r*  listed  in  Ubtf  U. 

Sup  1.  'n*  tytt*H«  or  ecuationt  it  to  vcd  using  me  full  tu(« 
*rc£uc**.  frew.  IK  «>:i«imiM  procedure.  Ihe  t/ttem- 
r#iv«v*  10  4  pitch  rate  iftpwlt*  of  1  deg/tec  It  Shewn 

i#s  f  i$  10. 

Step  2.  A’*.  f rc'v.'tO 4  to  live  direct  lift  control  t*rfact  *rt  set 
.o  /rro  to  del*  ratine  if  it*  direct  lift  control  la*  <4n  be 
#ll<Mft* i#4.  ‘he  resulting  AJCS  It  unstable.  Since  U.C 
It  *  pnrVpeCl*lcd  Control  and  IV  heavily  In  thf 

«S»V*  r#»f  anulyvlt.  U*  pr4CHC4l  t/tien  Mflll  *4»»  lo 
include  tilt. 

Slop  3.  Ik  tlneur  control  gains  were  deleted  to  Ctumine  If 
t#.e  «V«ator  cojld  to  elmtnauc.  1*4  numerical  results 
U  U&l*  H  indicate  that  elevato r  conirol  If  n&t  needed, 
hpwrver,  these  results  a*e  taved  on  rat  «alwet  and  C4A  bo 
misleading.  Ine  control  syttw  without  UK  elevator  ft 
«c  1**11/  slightly  Verdes s*d.  ‘he  system  response  to  4 
pitch  rau  i«tMh*  of  1  c*g/s«  It  iacmi  in  fig.  II.  Ire 
ef'ect  of  elector  control  it  guile  evident  when  fist.  10 
and  11  4 re  compared.  ?**  elevator  control  law  it  loft  in 
the  1/tlM. 

Sup  i.  J<ro  me  airspeed  feeeuatl  gam. 

Step  S.  Zero  me  potmen  feedPaca  saint  on  bom  control t. 

Stop  6.  Zero  the  wino  gust  intent!  1/  feeduc*  gains.  Th<*  system 

roprevontt  the  final  system,  before  adding  noise  y,  to  the 
glidetlope  deviation. 

Stop  7.  As  in  previous  cate*. ,  me  glideslcpt  deviation  it  no**  com¬ 
bined  with  the  fluctuation  nolte  and  men  filtered.  The 
filter  it  the**  in  ris.  0.  The  swe  tiro  cortuni  it  used 
at  in  the  patt  cate  **4re  If*.J  tec. 

Tne  final  system  save  very  good  results.  f-o  attempt  was  node  to 
improve  perforated*  through  addins  ;*‘ns  to  me  glides’©**  filter  or  by 
relating  rat  actuator  constraints  in  me  opiini/ation  routine.  In  an 
actual  application,  the  actuator  ratal  would  be  fined  and  the  design  would 
then  be  *ade  to  eoet  me  ratot  while  here  a  hypothetical  case  was  considered 
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Fig.  10.  Pitch  Rate  Impulse  Response  of  DLC  System 
With  Full  State  Feedback 
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Mg.  11.  Pitch  Rate  Impulse  Response  of  DIC  System 
Without  Elevator  Control 
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to  point  out  the  ease  with  which  the  optimal  approach  could  handle  this 
multi-input,  multi-output  s.ituation--not  to  achieve  the  best  possible 
performance. 

The  actual  numerical  values  for  the  feedback  gains  are  listed  in 
Table  XIV  for  the  case  with  airspeed  feedback,  and  in  Table  XV  for  the 
case  without  airspeed  feedback.  Both  tables  are  located  in  Appendix  A. 
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VII.  Summary  of  Results,  Conclusions,  and  Recommendations 


The  summary  of  results,  conclusions,  and  recommendations  reached 
in  determining  a  design  procedure  using  the  "optimal"  feedback  gains 
of  the  data  rate  analysis  are  enumerated  in  this  chapter. 

Summary  of  Results 

The  results  of  this  study  are  summarized  in  the  form  of  a  table 
and  five  time  histories  showing  aircraft  responses.  In  Table  X,  the 
PMA,  rms  elevator  actuator  rates,  rms  pitch  attitude,  and  rms  direct 
lift  control  rates  are  listed  for  the  "standard"  and  the  final  form  of 
each  of  the  four  cases  considered.  The  time  histories  (Figs.  12,  13, 

14,  15,  and  16)  show  the  system  response  to  an  initial  pitch  rate  impulse 
of  1  deg/sec  for  the  "standard"  and  for  each  case  considered.  The 
"standard"  and  the  two  cases  with  airspeed  have  approximately  the  same 
closed-loop  frequency.  The  two  cases  without  airspeed  seem  to  have  a 
higher  closed- loop  frequency. 

Conclusions 

1.  It  has  been  shown  that  a  set  of  optimal  feedback  gains  can  be 
the  basis  for  a  practical  control  system,  provided  the  following  considera¬ 
tions  are  properly  accounted  for  in  the  optimization: 

a.  The  cost  functional  reflects  system  performance. 

b.  Control  activity  is  realistically  accounted  for. 

c.  The  measurements  available  for  feedback  are  taken  into 
account  in  the  optimal  design. 
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TABLE  X 

Comparison  of  the  Final  System  for  Each  Case  and  the  Standard 
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♦This  value  was  computed  by  the  digital  computer  program  and  was  not  given  in  Ref  2. 
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Fig.  12.  Pitch  Rate  Impulse  Response  of  the  Standard  System 
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Fig.  13.  Pitch  Rate  Impulse  Response  of  the  DC-8 
With  Airspeed  Feedback 

AO 
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Fig.  <4.  Pitch  Rate  Impulse  Response  of  the  DC-8 
Without  Airspeed  Feedback 
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Fig.  15.  Pitch  Rate  Impulse  Response  of  the  Direct  Lift  Control 
DC-8  with  Airspeed  Feedback 
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Fig.  16.  Pitch  Rate  Impulse  Response  of  the  Direct  Lift  Control 
DC-8  without  Airspeed  Feedback 
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2.  A  systematic  procedure  was  developed  which  gives  a  good  "first 
cut"  at  a  practical  control  system. 

3.  The  procedure  is  "clean"  and  requires  no  "games”  or  "fiddling" 
with  the  feedback  gains. 

4.  The  AFCS  achieved  in  the  case  of  the  DC-8  without  airspeed 
feedback  compared  quite  well  with  the  standard  as  shown  in  Table  X  and  by 
comparing  Figs.  12  and  14.  If  the  high-frequency  oscillations  in  pitch 
attitude,  shown  in  Fig.  14,  are  considered  undesirable,  a  washout  on 
pitch  attitude  feedback  can  be  used  to  eliminate  th«r  (Ref  7). 

5.  The  incorporation  of  airspeed  feedback  in  the  control  law 
results  In  improved  approach  performance  (i.e.,  lower  PKA)  and  a  reduced 
sensitivity  to  changes  1r.  the  other  feedbacks.  When  interpreting  the 
results  in  Table  X,  the  reader  must  consider  the  value  of  the  ms  actuator 
rate  as  well  as  PHA  to  avoid  being  misled. 

6.  Dlnct  lift  control  was  handled  easily.  The  results  for  the 
two  DIC  cases  are  presented  in  Tabic  *  and  Figs.  15  and  16.  The  merit 

of  OIC  In  the  landing  approach  task  is  evident  from  the  results  In  Table  X. 

7.  The  rms  measures  of  "goodness"  are  good  for  a  "first  cut” 
design  but  are  not  sufficient  to  assure  a  desirable,  safe,  flyeble  AFCS. 
Further  analysis  il.e.,  pitch  rate  response  plots,  stability  tests,  phase 
and  gain  margin  checks,  etc.)  Is  required  to  assure  a  safe  system. 

Recommendations 

The  following  recommendations  are  suggested: 

1.  A  better  or  Improved  measure  of  landing  approach  performance 
is  needed.  PHA  is  not  adequate  even  when  rms  pitch  attitude  and  rms 
control  activity  are  considered.  What  is  needed  is  a  measure  of  ride 
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qualities  or  pilot  rating.  Riot  quality  can  be  '/value  ted  by  looking  at 
the  acceleration  spectra  at  the  cockpit  In  t’4  3  to  S  cps  region. 

2.  Adopt  the  approach  to  usks  other  than  landing  approach,  for 
Instance,  weapon  delivery. 

3.  Apply  this  procedure  to  design  an  automatic  lendirg  control 
t,ystn  for  raotely  piloted  vehicles. 
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W(kZ:*  A 

♦mental  Values  of  feedback  Gains 

ln«i  ifi/onMUOf.  in  this  apptndi*  is  included  for  the  bcntfi t  of 
anyone  wno  night  want  to  duplicate  this  study  or  do  further  research  In 
this  area.  All  values  presented  here  arc  Identical  to  the  values  used 
in  this  study.  The  feedbacks  retained  in  the  practical  system  are  marked 
v.1  th  an  asterisk  {*}. 

The  feedbacks  on  glideslopc  deviation  ar.d  airspeed  are  noved 
frrr*  Ve  pure  slates  (d,  Jd,  uAS)  in  the  full  state  feedbacc  system  to 
the  filtered  or  sensed  states  (df,  fd f,  us)  in  the  practical  system. 


67 


GGC/MA/73 -3 


TABLE  XI 

Feedback  Gains  for  the  Standard  System 
(data  rate,  continuous) 


State 

1 

1  Feedback  to  Elevator 

uas 

1  '  - 

1 

0.0 

w 

-  2.5576  xlO*2 

0 

5.8313 

q 

2.0000 

us 

1 

0.0 

w9 

0.0 

Of  1 

1 

-  1.0000 

yj 

0.0 

d 

0.0 

«e 

0.0 

/<* 

7.6803  *  10-* 

Of 

8.6700*  10-* 
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TABLE  XII 

feedback  Gains  for  the  Sub-optimal  System  with  Airspeed  Feedback 

(data  rate,  6  samples/sec) 


State 

Feedback  to  Elevator 

Feedback  to  Throttle 

uas 

3.3678  x  10's* 

-  .22531 

w 

-6.4191  x  10'3* 

9.5716  x 10-2 

e 

2.2657* 

-32.154 

q 

.85460* 

-13.897 

u9 

1.0442  xio*5 

-  .18910 

Wg 

-  3.2762  x  10  ** 

9.5802  xlO'3 

y» 

0.0 

1 

0.0 

d 

3.0483  x  TO"** 

-5.0314x10-* 

ae 

-  4.7493  x  10"2 

.78816 

J- 

2.6690  x  10-'** 

-2.3104x10-* 

6th 

2.1731  «10-' 

1 

-  2.3412  x 1 0“* 

i 


•This  feedback  Is  retained  In  the  practical  system. 
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TABLE  XIII 

Feedback  Gains  for  the  Sub-optimal  System  Without  Airspeed  Feedback 

(data  rate,  6  samples/sec) 


State 

Feedback  to  Elevator 

Feedback  to  Throttle 

uas 

4.8803  x  1  O'3 

-  1.3055  x  10-" 

w 

-  1.1965  x  10~2* 

6.2504  xlO’5 

e 

4.3265* 

-  2.1616  x 10-2 

q 

1.4094* 

-7.9522  x  10-* 

u9 

1  .2629  x  10‘3 

-1.1141  xlO“" 

Wg 

-  7.8283  xlO-- 

6.9985  xl0“6 

y2 

0.0 

0.0 

d 

7.4510  x  10"3* 

-4.3170x10-" 

«e 

-  7.6428  xlO"2 

4.4071  xlO*" 

J- 

3.4793  xlO“"* 

-  1 .3305  x  10  5 

6th 

5.7354  x  1 0-1* 

-2.4500  xlQ-‘ 

•This  feedback  Is  retained  In  the  practical  system. 
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TABLE  XIV 

Feedback  Gains  for  Sub-optimal  System  Without  Airspeed  Feedback 
Relaxed  Control  Activity  Constraints 
(data  rate.  6  samples/sec) 


State 

Feedback  to  Elevator 
(Gains  for  Practical  System  Only) 

uas 

0.0 

w 

-  2.1154  x lO"2 

e 

7.7203 

q 

2.1266 

U9 

0.0 

wg 

0.0 

! 

yz 

0.0 

d 

1.6108*10-2 

0.0 

fi 

4.6272  x lO"* 

Hh 

0.0 
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TABlE  XV 

Feedback  Gains  for  Direct  Lift  Control  System  with  Airspeed  Feedback 


State 

Feedback  to  Elevator 

Feedback  to  DLC  Surface 

uas 

6.4943  x 10“** 

-4.0792x10-** 

w 

-9.8753  x  10-"* 

7.1767  x10-** 

e 

.34236* 

-2.5057* 

q 

.14456* 

-  .98193* 

u9 

2.3002  x10-" 

-1.3089x10-* 

Wg 

-4.0206  x  10-* 

3.4499  x10“" 

yi 

0.0 

0.0 

d 

3.3161  x 10-"* 

-3.4505 xlO'3* 

6e 

-8.1636x  10-* 

5.4893  xlO-2 

fidlc 

3.2540  xlO“* 

-2.3556  x10-* 

*Th1s  feedback  Is  retained  In  the  practical  system. 
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TABLE  XVI 


Feedback  Gains  for  Direct  Lift  Control  System  Without  Airspeed  Feedback 


State 

Feedback  to  Elevator 

Feedback  to  DLC  Surface 

lJas 

1.8430  xlO'4 

-5.4850  x  10-* 

w 

-3.5975  x  10“*** 

1.4798  x10-** 

e 

-  .12807* 

-5.2815* 

q 

5.2266  xlO'2* 

-1.7123* 

ug 

4.9343  x  10"* 

-  .10923 

w9 

-1.3738 x 10-s 

9.5572  x  10-- 

y? 

0.0 

0.0 

d 

I 

1.3172  xlO--* 

-9.3767  x 10--* 

ie 

-2.9351  xlO'* 

5.2860  xlO'2 

6dlc 

1.1819x10-* 

-4.8113  xl0“2 

♦This  feedback  Is  retained  In  the  practical  system. 
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